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We present phase-resolved spectroscopy and photometry of a source discovered 
with the Chandra Galactic Bulge Survey (GBS), CXOGBS J174444.7- 260330 (aka 
CX93 and CX153 in the previously published GBS list). We find two possible values 
for the orbital period P, differing from each other by ~13 seconds. The most likely 
■^j- solution is P = 5.690f 4(6) hours. The optical lightcurves show ellipsoidal modulations, 

whose modeling provides an inclination of 32±1° for the most likely P. The spectra are 
t-H dominated by a K5 Vcompanion star (the disc veiling is <5%). Broad and structured 

emission from the Balmer lines is also detected, as well as fainter emission from Hel. 
From the absorption lines we measure K2 =117±8km/s and wsini = 69 ± 7km/s. By 
I solving the system mass function we find Mi =0.8±0.2M Q for the favored P and i, 

£NJ consistent with a white dwarf accretor, and M2 =0.6±0.2 M . We estimate a distance 

t— I in the range 400 — 700 pc. Although in a low accretion state, both spectroscopy and 

photometry provide evidence of variability on a timescale of months or faster. Besides 
. *-H finding a new, long orbital period cataclysmic variable in a low accretion state, this 

work shows that the design of the GBS works efficiently to find accreting X-ray binaries 
5— 1 in quiescence, highlighting that the spectra of CVs in a low-accretion state can at times 

appear suggestive of a quiescent neutron star or a black hole system. 

Key words: stars: individual (CXOGBS J174444.7-260330) — individual (CX93 
=CXf 53) — accretion: accretion discs — X-rays: binaries — (stars:) novae, cataclysmic 
variables — stars: dwarf novae 



1 INTRODUCTION 

We present optical follow-up of the X-ray source CX- 
OGBS J174444.7-260330 J _d^overed in the Galactic Bulge 
survey (GBS) of |Jonker et"aLlpoTl) . For the ease of read- 



ability, we will refer to the source with its label in the GBS 
source list, CX93. One of the main goals of the GBS is to 
detect quiescent X-ray binaries (XRBs) - namely binary sys- 



* data from ESO programs 085.D-0441(A) and 087.D-0596(D) 
f email : e.m.ratti@sron.nl 



tems where a white dwarf (WD), a neutron star (NS) or 
a black hole (BH) is accreting matter from a companion 
star - that are suitable for dynamical measurements of the 
mass of the compact object. Those mass measurements can 
be done through phase-resolved optical spectroscopy of the 
donor star, provided that good constraints can be put on 
the binary inclination. If the system inclination i is known, 
the mass of the accreting compact object can be obtained by 
solving the system mass function (see, for instance, |Charles| 
fc Coe|2006[ ). This requires measuring the orbital period P, 



the semi-amplitude of the radial velocity curve of the com- 
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Figure 1. 20"x20" finding charts showing the near infrared (J- 
band) and the optical (R-band) counterpart to CX93, indicated 
by the black arrows. The left panel is from the variable sources in 
the Via Lactea (VVV) survey, where the VVV source detections 
are indicated by circles. Two overlapping circles at the counter- 
part position correspond to the Chandra position and to that of 
the NIR source (R.A. =17 h 44 m 44!790, Dec. = -26° 03' 30'.'08, 
with an accuracy better than 0.1"). The right panel is an image 
from the Visible Multi-Object Spectrograph (VIMOS). 



panion star K2 and the ratio between the mass of the donor 
and that of the accretor, q = M2/M1. In Roche lobe filling 
XRBs, with tidally locked companion stars, q is a function of 
K2 and of the projected rotational velocity of the secondary 



star vs'mi (Wade & Home 19881, which is measured from 



the broadening of the stellar absorption lines. 
Prime targets for this study are eclipsing systems, where the 
inclination can be derived from the eclipse duration based 
on geometrical arguments only. Among non-eclipsing sys- 
tems, the inclination can be estimated by modeling the ellip- 
soidal variation in the optical lightcurve, caused by the dis- 
tortion of the companion star shape associated with Roche- 
lobe overflow (care must be taken if other continuum sources 
contribute to the lightcurves). 

Because CXOGBS J174444.7-260330 was initially not 
identified as a duplicate, it has two labels in the GBS source 
list, CX93 and CX153 (the offset between the two detections 
is 3"). The best X-ray source position, from observation 
CX93, is R.A. =17 h 44 m 44!791(9), Dec. = -26° 03' 3(/.'3(l) 
(the la uncertainties, indicated in between brackets, do not 
include the 90% confidence 0'.'6 Chandra boresight error). 
The source was identified as a quiescent X-ray binary can- 
didate based on preliminary low-resolution spectra of the 
optical counterpart (Figure [T]). A near infrared (NIR) coun- 
terpart to CX93 was found in the data from the VISTA 
(Visible and Infrared Survey Telescope for Astronomy) vari- 
able sources in the Via Lactea (VVV) survey (Minniti et al. 
2009, Catelan et al. 2011, Saito et al 2011), with magni- 
tude 14.81±0.02 in the J-band, 14.01±0.02 in the H-band 
and 13.76±0.03 in the K-band. Here we present phase- 
resolved photometry and optical spectroscopy of the target 
performed with several instruments. We detected signatures 
from the mass donor star, which allowed us to measure the 
orbital ephemeris and constrain the mass of the accreting 
compact object. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Phase-resolved imaging 

CX93 was observed in 2010 Jul. using the Mosaic-2 in- 
strument on the Cerro Tololo Interamerican Observatory 
4 m Blanco telescope (Table [TJ . Thirty three 120s-long ex- 
posures were collected through the Sloan r' filter. Initial 
data processing was done by the National Optical Astron- 
omy Observatory (NOAO) Mosaic pipeline ( |Shaw|[2"009| ). 
The lightcurves were extracted using the ISIS image sub- 
traction code flAlard fc Lupton|[l998l |Alard||1999|, which 



provided clean subtracted images. A variable optical coun- 
terpart coincident with the X-ray position of CX93 was 
clearly identified, despite the fact the variability is only at 
the level of ~5%. We independently analyzed the images 
using DAOPhot II (Stetson 1987) and obtained compara- 
ble lightcurves, although of somewhat poorer quality than 
those obtained with ISIS. Since ISIS only yields differential 
count rates, we used DAOPhot II on the best reference im- 
age to measure the baseline (non-variable) count rate. This 
together with the differential count rates from ISIS yielded 
total count rates as a function of time. We finally converted 
these to approximate magnitudes using conversion factors 
supplied by the NOAO pipeline. These are deduced by com- 
parison to USNO B1.0 stars, and have an estimated uncer- 
tainty of at least 0.5 magnitudes. The magnitude of CX93 is 
17±0.5 magnitudes in the r'-band. 

CX93 was also observed using the Direct CCD Camera on 
the Henrietta Swope Telescope at Las Campanas Observa- 
tory on the night of 2011 Jun. 27/28. A sequence of 60 sec 
exposures was obtained through a Gunn r filter for about 
six hours. Bias correction and flat fielding was performed 
through standard CCD data processing in iRArQ ISIS was 
used for image subtraction photometry as for the Mosaic-2 
data. 



2.2 Optical spectroscopy 

We combined spectroscopic observations from a number of 
telescopes: the Visible Multi-Object Spectrograph (VIMOS) 
and the Focal Reducer and Spectrograph (FORS) at the 
Very Large Telescope, the Inamori Magellan Areal Cam- 
era and Spectrograph (IMACS) at the Magellan telescope, 
the ESO Faint Object Spectrograph and Camera (EFOSC) 
at the New Technology Telescope (NTT) and the Magel- 
lan Echelle (MagE) spectrograph, also at the Magellan tele- 
scope. A list of the observations and the instrument set- 
tings is given in Table [T] The NTT/EFOSC data are not 
included in the table and in the rest of this work because of 
the poor quality of the data due to weather conditions. Tem- 
plate spectra of the stars HD163197 (spectral type K4IV) 
and HD130992 (K3.5V) were also observed with FORS and 
MagE, respectively. The images were corrected for bias, 
flat-fielded and extracted using the Figaro package within 
the Starlink software suite and the packages Pamela and 
Molly developed by T. Marsh. We used sky flats for the 



1 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities 
for Research in Astronomy (AURA) under cooperative agreement 
with the National Science Foundation. 
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Table 1. Journal of the observations. The resolution of the spectra is measured from the width of the night sky lines. 



Instrument 


Date 


# 


Exposure 


Grism/Filter 


Seeing 


Slit width 


Binning 


Sp. range 


Resolution 


Dispersion 




CUT) 




to 




(arcsec) 


(arcsec) 




(A) 


(A) 


(A/px) 


Spectroscopy 


VIMOS 


2011 Apr. 29 


2 


875 


MR-2.2 


0.6 


1.0 


lxl 


4879-10018 


10 


2.50 


IMACS 


2011 Jun. 23-26 


i 


1200-600 


Gri-300-17.5 


0.8-1.2 


1.0 


1x1 


4365-6654 


4.5 


1.28 


FORS 


2011 Jul. 4 


i 


900 


1200R+93 


0.9-1.4 


1.0 


2x2 


5957-7280 


2 


0.38 


MagE 


2011 Aug. 6 


1 


900 


Echelle 


1.2 


1.0 


lxl 


3100-11200 


1.6 


0.47 


Photometry 


Mosaic 


2010 Jul. 12-18 


33 


120 


Sloan r' 


1.3 












Swope 


2011 Jun. 27-28 


360 


60 


Gunn r 


1.5 













flat-fielding and we subtracted the sky background by fit- 
ting clean sky regions along the slit with a second order 
polynomial. The spectra were optimally extracted follow- 



ing the algorithm of Home ( 1986 ) implemented in Pamela 



and wavelength-calibrated in Molly with a final accuracy 
of 0.02 A for FORS (using arc exposures taken during day- 
time), 0.04 A for IMACS, 0.09 A for VIMOS and ~0.04 A 
for MagE. In each spectrum, the wavelength calibration was 
checked and corrected for shifts with respect to the position 
of the night-sky OI lines at 5577.338 A and/or at 6300.304 A 



(Osterbrock et al. 1996). Each spectrum has been normalised 



dividing by a spline fit to the continuum, with a maximum 
order of 10. The full MagE spectrum was extracted, but only 
order 9 (5630-6673 A) and 10 (6265-7415 A) were used for 
our dynamical study of CX93, as they cover the spectral re- 
gion of the Ha and Cal, which is rich in absorption features 
for K type stars and where the signal-to-noise ratio in the 
spectrum is highest. The MagE spectra show a complex con- 
tinuum, mainly caused by instrumental response variations 
across the relevant orders. In order to normalise the spec- 
tra, we clipped two ~300 A-long pieces from the order 9 
spectrum (5920-6240 A and 6310-6630 A), and one from the 
order 10 (6300-6750 A). The regions were selected to keep 
as many absorption lines as possible in one spectrum. By 
selecting short pieces from the spectrum we could achieve 
a good fit to the continuum on each piece with a relatively 
low-order spline (order 10), avoiding as much as possible 
the risk of altering the equivalent width of the lines that 
occurs when using high order splines. The overlap between 
the orders 9 and 10 allows for a double-check of the results. 



A is next to a blend of the Nal doublet lines in absorp- 
tion (A5895.92, A5889.95) mainly produced by the interstel- 
lar medium. The Mgl triplet at A5167, A5172 and A5183 
is visible, although not resolved. The Cal triplet (A6102, 
A6122 and A6162) and a blend of stellar lines in absorption 
at A6495 are also present. 

The Doppler-corrected average of the FORS spectra, in the 
central panel of Figure[2] shows the region of Ha, with better 
resolution than that of the IMACS spectrum. Hel in emis- 
sion at A6678.149 and a forest of absorption lines are visible, 
in particular from Fell at A6432.65 and A6457 and Fel at 
A6430.85. Finally, the bottom panel of Figure [2] shows a sec- 
tion of the Doppler-corrected average of the VIMOS spectra, 
displaying the main features in the red part of the spectrum 
of CX93: the Call triplet in absorption, at A8498, A8542 and 
A8662, and a number of emission lines from the Paschen se- 
ries of hydrogen. A telluric molecular band of water is also 
present. 

In most of the spectra, the Ha emission line appears single- 
peaked but composite, with a main peak overlapping with 
at least one side wing. Both components are broad, with 
FWHM of the order of 350-450 km s" 1 . The equivalent width 
(EW) of the line is roughly constant around ~ —13/16 A 
across the IMACS and FORS spectra, with no dependence 
on phase. In the MagE spectrum, taken one month later, it 
drops significantly to -2.9±0.7 A. The two VIMOS spectra, 
collected two months earlier than all the other data-sets, 
suggest a more active state of the source compared to the 
IMACS and FORS observations. The EW is -133±1 A and 
-140±1 A in the two spectra. The H/3 line has similar EW 
to Ha and an even larger FWHM of ~ 700 km s -1 . 



3 ANALYSIS AND RESULTS 



3.2 Binary period 



3.1 Spectral features 

The spectra of CX93 display a number of absorption lines 
from the stellar atmosphere of the companion star, but also 
strong emission lines consistent with Ha and H/3 and with 
Hel at A6678.149 and A5875.618. For each instrument, we 
combined all the spectra into an average spectrum, after 
correcting for the Doppler shift of the lines due to the or- 



bital motion (see Section 3.3 for the measurement of the 
orbital shifts). The top panel of Figure [2] shows the average 
of the IMACS spectra where the Ha and H/3 lines appear 
as prominent emission features. Hel in emission at 5875.618 



The orbital period of CX93 was determined using the 
Mosaic-2 lightcurve, combined with spectroscopic informa- 
tion. The Swope data provide only loose constraints on the 
system parameters, which we therefore only used as a con- 
sistency check of our results. In order to obtain the best 
possible accuracy on P, we followed three main steps. 

• An initial Lomb-Scargle periodogram of the Mosaic-2 
data revealed an apparent period of 0.12 d. Some disper- 
sion was seen when the Mosaic-2 lightcurve was folded on 
this period indicating that this is the first harmonic. Folding 
on twice this period yielded an asymmetric double-humped 
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Figure 2. Main spectral features from CX93 in the IMACS 
(top panel), FORS (middle panel) and VIMOS (bottom panel) 
Doppler-corrected average spectra. 



lightcurve (Figure[3] top panel). A period of P=0. 23710(5) d 
was measured by fitting the lightcurve with a sum of two 
sine waves, with a 2:1 frequency ratio, and amplitudes and 
relative phases allowed to vary freely. Error bars were esti- 
mated in the usual way from \ 2 totals, after adding an ad- 
ditional error to the Mosaic-2 datapoints to represent unre- 
solved flickering. The Swope data provide a less constrained 
P =0.23(1) d, consistent with the Mosaic-2 result. Unfortu- 
nately a joint fit to both data sets was not well constrained 
due to the difficulty in extrapolating the ephemerides across 




0.5 1.0 1.5 

orbital phase (two cycles shown) 




Figure 3. Top panel: Mosaic-2 lightcurve. The solid line repre- 
sents the best fit model assuming ellipsoidal variations plus a 
single bright spot. Bottom panel: LCO lightcurve, fitted with 
a pure ellipsoidal modulation model for an inclination of i = 
32.3°(dashed line) and i = 36.2°(drawn line). The observations 
are binned for clarity, but fitting was done on unbinncd data. 



a year (see below), and possible changes in lightcurve mor- 
phology. 

• The Mosaic-2 P was then used to phase fold the spec- 
troscopic observations and construct an initial radial veloc- 
ity curve (rvc, the procedure is described in Section [3. 3| for 
the final orbital period). The rvc looked consistent with or- 
bital motion of the companion star, and was fitted with a 
sine wave, measuring To in between the time of the spec- 
troscopic observations. We double-checked the measurement 
of the orbital period leaving it free in the fit, obtaining 
P = 0.23710(5) d as expected. 

• At last, To was compared with the time of the phase 
in the Mosaic-2 lightcurve, Tm. We refined the measure- 
ment of P by considering that it must be such that the 
difference To— Taj corresponds to an integer number of or- 
bital cycles. Both the minima in the Mosaic-2 lightcurve po- 
tentially correspond to Tm, since the P measured above is 
not accurate enough to unambiguously phase the lightcurve, 
collected one year earlier than the spectroscopic observa- 
tions. Depending on the phasing, we find three possible so- 
lutions for P. If Tm is at the deepest, primary minimum 
of the Mosaic-2 lightcurve, then P = 0.237089(3) d, with 
To — Tm = 14898 x P. If Tm corresponds to the secondary 
minimum, than P =0.237169(3) d with T -T A f = 14888 xP 
or P = 0.237009(3) d with T -T A f = 14889 x P. The rvc 
fitting favors the first solution, as it is the closest to the 
best-fitting period. The second solution is less likely but, 
yielding a A\ 2 of 3.5 with respect to the best-fit x 2 > it has a 
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probability of ~ 20% to be correct. The last solution is very 
unlikely, yielding A^ 2 ~ 15. 

In conclusion, P = 0.237089(3) d is the statistically most 
likely orbital period for CX93, although P = 0.237169(3) d 
can not be ruled out at high confidence. 



3.3 Radial velocity curve 

The orbital Doppler shifts of the companion star in CX93 
were measured from the absorption lines in the spectra, by 
cross-correlating with the K4 IV template star spectrum ac- 
quired with FORS. The emission lines were masked. Since 
we were using observations from different instruments, the 
template and each target spectrum were re-binned to the 
same dispersion and broadened to the same spectral reso- 
lution before the cross-correlation. For MagE, the order 10 
and the two parts in which the order 9 was split were cross- 
correlated separately with the FORS template, and the re- 
sults were averaged in order to reduce the uncertainty. For 
all the spectra, the uncertainty on the velocity offsets due 
to the wavelength calibration has been added in quadra- 
ture to the velocity errors from the cross-correlation. The 
observations were phase folded with the most likely orbital 
period from Section |3.2| obtaining the rvc shown Figure [4] 
(top panel). 

The rvc was fitted with a circular orbit of the form 
v(^>)=7+K2 sin(27r0 + ip), providing a large reduced \ 2 01 
~ 13 (8 degrees of freedom). The uncertainties on the pa- 
rameters were estimated assuming that the sinusoidal model 
was correct, and multiplying the errors of the individual ve- 
locity shifts by 3 to reach a reduced x 2 close to 1. To was 
calculated near the middle of the time of the observations 
and so that <p = 0, according to the convention of phase at 
inferior conjunction of the companion star. The period was 
fixed in the fit. The resulting parameter values are: 

T = 2455743.1892 ± 0.0029 HJD/UTC 

K 2 = 117 ± 6 km s" 1 
and 7 = —69 ± 6kms _1 , in the rest frame of template star 
employed for the cross correlation. A systemic radial veloc- 
ity of 29.8±0.4kms _1 was found for the latter, by fitting a 
Gaussian function to the Ha absorption line and measuring 
the offset of the line centroid with respect to its rest-frame 
wavelength. The systemic velocity for CX93 is thus 

7 = -39±6 kms" 1 . 
Figure |4]( bottom panel) shows the rvc for the Ha line, where 
the orbital velocities were measured through Gaussian fit- 
ting to the line peak in each spectrum. The curve was con- 
structed and fitted as we did above. The orbital motion is 
close to that of the companion star (from the absorption 
lines), only slightly leading in phase by tp — 0.08±0.03. The 
fit provides K em = 58 ± 6kms~ 1 , smaller than the radial 
velocity semi-amplitude we measured from the absorption 
lines, suggesting an origin of the Ha line towards the in- 
ner face of the companion. The CX93 systemic velocity is 
—45 ± 4 kms" 1 , consistent at the la level with the result 
above. Similarly, we measured radial velocity shifts for the 
Ha wing feature. We find variable velocity shifts from ~ 10 
up to ~600 km s _1 towards the red or the blue, with no clear 
dependence on the orbital phase. 




0.4 0.6 
Orbital phase 



Figure 4. Radial velocity curve measured from the absorption 
lines (top panel) and from the peak component of the Ha emission 
line (bottom panel). The data are fitted with a circular orbit (solid 
line), increasing the error bars in order to obtain a reduced \ 2 
close to 1. 
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Figure 5. Optimal subtraction of the spectrum of the stellar 
template HD13099 from the spectrum of CX93, both acquired 
with MagE, providing the measurement of vsini. The image 
shows the region around Ha from MagE order 9. The emission 
lines were masked in the optimal subtraction (see text Section 

E21. 



3.4 Rotational broadening and mass ratio 

The observed FWHM of the absorption lines in the spec- 
tra is determined by the intrinsic line width, expected to 
be dominated by rotational broadening, broadened by the 
instrumental resolution profile and smeared by the motion 
of the companion star during the integration time of one 
observation. In order to measure the intrinsic line width, 
namely usini, we analysed the MagE spectra, which have 
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the best resolution among our data-sets. We considered all 
the absorption lines in the spectrum (none of which is known 
to be strongly affected by thermal or pressure broadening), 
masking the emission lines. The spectral type of the MagE 
template, K3.5V, is likely close to that of the companion star 



in CX93 (see Section |3J5j). Steegh s fc Jonker| ( |2007| > showed 
that a small difference in the spectral type of target and 
template star do not significantly effect the measurement of 
v sin i in close binaries. The integration time smearing on the 
target spectra was taken into account by artificially smearing 
the template spectrum by 27rTif cos(27r</!>)/Pkms -1 , where 
T is the duration of one exposure on CX93. After that, we 
broadened the template spectrum with different trial values 
of -usini. For each trial value, the template was subtracted 
from the target spectrum and a \ 2 test was performed on 
the residuals. Different scalings of the template are tried 
to account for possible disc veiling The broadening that 
gives the minimum x 2 is a measure of the actual vsini. A 
value of 0.5 for the limb darkening was assumed, appropriate 
for a main sequence star as we expect for a ~5 hour orbit. 



Steeghs & Jonker ( 2007 1 tested the effect of the uncertainty 



on the limb darkening coefficient on vs'mi, showing that it 
is ~ f kms~ 1 for a difference of 0.25 in the limb darkening. 
This is negligible with respect to the error on our measure- 
ment (see below). Finally, since we subtract template and 
target spectra acquired with the same instrument, the in- 
strumental resolution profile is not affecting the procedure. 
The uncertainty on v sin i was estimated following the Monte 



Carlo approach used by Steeghs & Jonker (20071. We copied 



the target spectrum 500 times, using a bootstrap technique 
where the input spectrum is resampled by randomly se- 
lecting data points from it. The bootstrapping maintains 
the total number of data points in the spectrum. For each 
bootstrap copy, one value of v sin i was measured, corre- 
sponding to the minimum \ 2 of the residuals as described 
above. The distribution of «sini obtained from the 500 
copies is well described by a Gaussian, whose mean and root- 
mean-square (rms) provides the best-fit ?;sini and its la 
error. Our final measurement of the rotational broadening 
is the uncertainty-weighted average of the measures from 
the first and second piece of the MagE order 9, usini = 
69 ± 7 km s" 1 (Figure |5| . 

Given v sin i and K2 we could calculate the system mass ra- 
tio q = Mi I Mi from the relation 



0.462<? 1/3 (l + g) 



\2/3 



( |Wade fc Horne|f988[ ) . The uncertainty on q was again esti- 
mated with a Monte-Carlo simulation, calculating q for f 000 
sets of v sin i and K2 randomly selected within the la uncer- 
tainty on the parameters. The average of the f 000 simulated 
measurements provides q = 0.7±0.f, where the uncertainty 
is the standard deviation of the sample. 



2 the optimal subtraction is performed using the OPTSUB com- 
mand in MOLLY. The routine works on normalized spectra, as- 
suming Fd = F a — Ft * /, where F s is the flux of the target, Ft 
of the template, F^ is the disc contribution outside the lines and 
< / < 1 is a constant, the optimum factor. The routine finds the 
/ that minimizes the structure in the residuals between F4 and a 
smooth version of F^ (< Fj, >), accounting for broad features in 
the disc spectrum. The disc veiling is 1 — /. 



3.5 Companion star spectral type 

In order to constrain the spectral type of the mass donor in 
CX93, we considered a set of 21 high signal-to-noise template 
spectra of G and K stars of luminosity class V, III and IV, 
acquired with the William Herschel and the Isaac Newton 
telescopes between 1992 and 1994. Some of them were pre- 
viously used for spectral type classification by |Casares et al.] 
\im&\ and |Torres et al.| ( |2002[ ). We subtracted each tem- 
plate spectrum from the Doppler-corrected average of our 
FORS observations, performing a % 2 test on the residuals of 
the subtraction. The template resulting in the minimum x 2 
provides our best estimate for the source spectral type. The 
spectra were shifted to the rest frame of the FORS average 
spectrum, and degraded to match the FORS sampling and 
the line broadening in the target spectra. The appropriate 
broadening (accounting for the difference in resolution be- 
tween the templates and the FORS spectrum and the line 
broadening intrinsic to the source) was found in a similar 
manner as we did for measuring usini: different broadened 
versions of each template were subtracted (through optimal 
subtraction, accounting for veiling) from the CX93 spec- 
trum, until the minimum \ 2 was obtained from the fitting 
of the residuals. As for usini, we assumed a limb darkening 
of 0.5. 

The minimum \ 2 is obtained for the K5 V template 
6ICYGA, with an optimum factor of 0.98±0.01, indicating 
a disc veiling of ~ 2%. It is reasonable to assume an uncer- 
tainty of one spectral type when using this procedure. 



Using the same method decribed in Ratti et al. (20101 (but 



adopting standard stellar magnitude in the 2mass J, H and 
K filters from Mamajek 201]Fj we find that the spectral type 
K5 V is fully consistent with the J, H and K magnitude of 
CX93 from VVV, for a reasonable value of the extinction of 
Ay = 2.0 ± 0.6 (Ax = 0.23 ± 0.07, see Discussion). 

3.6 Ellipsoidal modulations and system inclination 

The phase folded Mosaic-2 and Swope lightcurves in Figure 
[3] show variability between the two data-sets. The Swope 
lightcurve displays pure ellipsoidal modulations, with equal 
maxima and the deepest minimum at phase 0.5, while the 
Mosaic-2 one, collected one year earlier, presents asymmet- 
ric maxima. 

We attempted to determine the system inclination by mod- 
eling the lightcurves with the XRbinary program written by 
E.L. RobinsorQ The Mosaic-2 Sloan r' and the Swope Gunn 
r filters were represented by a 5548-6952 A and by a 6181- 
7212 A square bandpass, respectively. In the modeling we 
assume K2 and q as determined in Sections |3.3| an d |3.4| and 
that the mass donor is a K5V star (see Section |3.5[ ) with 
an effective temperature of 4500 K (Mamajek 2011). As we 
did in Section [3. 2| the Mosaic-2 data were primely used to 
constraint the inclination, while the Swope observations was 
kept for consistency checks only, due to the poorer quality 
of the data. 

Depending on the phasing, we fitted the Mosaic-2 lightcurve 
(Figure [3] top panel) with two models. For the most likely 

3 http://www.pas.rochester.edu/~emamajek/ 
EEM_dwarLUBVIJHK_colors_Teff.dat 

4 http://pisces.as.utexas.edu/robinson/XRbinary.pdf 
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period, the lightcurve has the deepest minimum at phase 
0.5 and the highest maximum at phase 0.25. We modeled 
this with ellipsoidal modulations plus a single (corotating) 
hot spot on the outer edge of an otherwise non-luminous 
disc, centred so that it faces earth at phase 0.25 (which is 
unusual, see Discussion). This provides a reduced \ 2 01 ~4 
(30 d.o.f.) and an inclination of 32.3 ± 1.3°. 
For the alternative orbital period, the minimum at phase 0.5 
is no longer the deepest, which implies heating of the side 
of the secondary facing the primary. We modeled this as- 
suming irradiation from the primary component. On top of 
ellipsoidal variations, a disc-edge spot, now at phase 0.75, is 
again included in the model to account for the asymmetric 
maxima. Since this model depends on the (unknown) disc 
geometry, two disc versions were tried, with a height to outer 
radius ratio of 1 to 12 and 1 to 24. The resulting system in- 
clination is 36.2 ± 1.3° and 34.5 ± 1.3°, respectively. 
The Swope data, modeled with ellipsoidal modulations only, 
are consistent with all the values of the inclination deter- 
mined above. Because the Swope observations were taken 
close to the spectroscopic ones (one night after the last 
IMACS pointing), and the disc veiling was ~2% in the spec- 
tra, it is safe to assume very little disc contribution to the 
Swope lightcurve. Its consistency with the inclination esti- 
mated from the Mosaic-2 data therefore suggests that the 
latter are not heavily affected by an accretion disc contri- 
bution either (on the effect of a disc contamination on i see 
Cantrell et aTpoTo i. 



3.7 Mass of the stellar components 

The mass function /(Ma) : M\- 



PK-2 



0.04 of 



(1+«F 2ttG 

CX93 is too small to constrain the nature of the primary. 
However, including our measurements of i and q, we can 
solve it for all dynamical parameters. For the most probable 
period (see Section 3.2 I, we obtain Mi = 0.8 ± 0.2 Mq and 
M2 = 0.6 ± 0.2 Mq. The uncertainty was obtained with a 



Monte-Carlo method, as we did for q in Section |3.4| For the 
alternative period the inclination is more uncertain and the 
masses are reduced to Mi ~ 0.6 Mq and M2 ~ 0.4 Mq. 



3.8 UV counterpart 

We searched for an ultraviolet (UV) counterpart to CX93 
in 900 seconds long observations from the Galaxy Evolution 
Explorer (GALEX). The closest object to the source is a 
2<r-detection at an angular distance of 2.2" from our Chan- 
dra position. Even assuming that the detection was real, the 
position is not consistent with that of CX93 (the ler accu- 
racy of the astrometry of GALEX is 0'.'4, plus a systematic 
uncertainty of 0'.'2 in Dec). The non-detection of a counter- 
part in GALEX provides an upper limit to the magnitude of 
CX93 of ~22 AB magnitudes in the near UV GALEX filter 
(at 2500 A, limiting magnitude at 3<r). 



4 DISCUSSION 

We have performed phase resolved optical spectroscopy and 
photometry of the optical counterpart to CX93, constrain- 
ing the system parameters. We found an orbital period of 



more than 5.6 hours (most likely P = 5.69014(6) hours) 
and an inclination of 32-36° (most likely 32.3 ± 1.3°). The 
best-estimated masses of the primary and secondary star 
are O.8±O.2M0 and O.6±O.2M0, although lower masses 
are possible depending on the phasing of the Mosaic-2 
lightcurve. The latter displays minima consistent with el- 
lipsoidal modulations, but asymmetric maxima (this is of- 
ten called O'Connel effect, see |Wilsey fe Beaky||2009l for 
a review). Two values of the period, differing by ~13 sec- 
onds, provide an acceptable fit to both the lightcurve and 
the radial velocity data, causing an ambiguity of half cycle 
in the phasing. For the most likely period, we could model 
the lightcurve with ellipsoidal modulations plus a disc spot 
accounting for the highest maximum, at phase 0.25. The 
phasing of the disc spot is unusual though, as a mass trans- 
fer stream hitting the accretion disc is expected to produce 
a hot spot that leads the donor. Alternative effects causing 
the shape of the lightcurve, consistent with the phasing, are 
tidal interaction of the companion star with the outer disc 



(Frank et al. 2002 although the most commonly observed 



tidal interaction effect, superhumps, can only happen in sys- 
tems with q < 0.25, |Hynes||2012|) or starspots reducin g the 
light at one maximum (see, e.g., Wilsey fc Beaky|2009 1 . The 
latter are often invoked to explain to the O'Connel effect, 
but a realistic model of starspots is hard to produce with 
lightcurve modeling codes, especially when a single band 
lightcurve is available. Compared to the disc-spot model, the 
effect of a starspot model on the system parameters would 
be to increase the inclination, leading to smaller masses. 
All our data indicate that the primary is a WD, placing 
CX93 in the class of cataclysmic variables (CVs). The mass 
of the donor is consistent with the spectral type favored 
by our spectroscopic observations (K5V, typical mass is 
O.67M0; [D"rilling fc Landolt|2000| , or slightly under -massive 
for the highest inclination scenario. Under-massive compan- 
ions are often observed in XRBs and CVs as a result of their 



accretion history (e.g., Her X-l and Cyg X-2, see Tauris & 
van den Heuvel|2006 Knigge 2006). The de-reddened upper 



limit to the UV flux that we could obtain from the Galex ob- 



Cardelli et al. 



Koester et al. 



servations (assuming Ay as derived in section 3.5 \ is rather 
uncertain, due to possible deviation of the extinction law 
compared to the standard one we used (from 
|1989[). Still, by comparing WD models from 
(2005) with our UV upper limit, we conclude that the non- 
detection of a UV counterpart is consistent with the pres- 
ence of a 15000-30000 K WD (in agreement with the WD 
temperature we find for the irradiation model of the Mosaic- 
2 lightcurve) while a much hotter object would have been 
detected. 



4.1 Distance and X— ray luminosity 

By comparing the typical K-band magnitude of a K5V 
star (M_r- =4.42, Mamajek 2011) with the observed mag- 
nitude rriK = 13.76 ± 0.03, and assuming Ay = 2.0 ± 0.6 
(Ak = 0.23 ± 0.07) as estimated in Section 3.5 we obtain 
the distance d =664±23pc to CX93. This estimate relies 
on the assumption that the absolute magnitude of the com- 
panion is typical of a normal main sequence star, which is 
not true for CVs in general (Knigge 2006). The typical ra- 
dius of a K5V star is 0.72 Kq prilling fc Landolt||2000|) 
Using the mass-to-radius relation found by Kniggel (2006) 
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for secondary stars in long orbital period CVs, we find that 
our best mass estimate for the donor, 0.6±0.2Mq, yields a 
radius of 0.62 ± 0.07 R© (this is consistent with the Roche 
lobe radius of the companion Rrl2 =0.62 Rq for Mi = 0.8 
and q = 0.7). The secondary in CX93 is then fainter than a 
typical K5V by 0.1 to 0.6 magnitude, allowing a distance of 
~500 pc. Further uncertainty to the distance comes from the 
assumed spectral type. Allowing the spectral type to be off 
by one, our observations suggest a likely distance between 
~ 400 - 700 pc. 

Although seemingly high, the extinction we estimated is con- 
sistent with a nearby object. In fact, extinction maps ob- 
tained from red clump stars within VVV, with a resolution 
of 2'x2', indicate E(B-V)=1.67 (Ay ~ 5.1) in the direction 
of CX93 ( [Gonzalez et al.|2011 1. 

Converting Ay into an hydrogen column density (following 



Giiver & Ozel 



20091 yields N H = (0.4 ±0.1) 



Assuming a power-law spectrum with index 1.6, the 14 
counts observed by Chandra from CX93 in the full ACIS 
range (0.3-8 keV) correspond to an unabsorbed flux of 
3.9xl0" 13 ergs -1 cm -2 . For d =664 pc the 0.3-8 keV X-ray 
luminosity of CX93 is 1.8x 10 31 ergs -1 . Compared with typ- 
ical CVs luminosities (~ 10 3 
with a low accretion state. 



- 10 34 erg s 1 ) this is consistent 



4.2 The difficult classification of quiescent XRBs 

The low X-ray luminosity, the the lack of clear signatures 
of ongoing accretion and the little disc veiling in the opti- 
cal spectra indicate that CX93 was in a low state during 
most of our observations. Disc emission is possibly respon- 
sible for the wings at the base of the Ha line, showing high 
radial velocities of a few hundred kms -1 , but the peak of 
the line traces the donor star (only slightly leading in phase) 
and can hardly be reconciled with disc motion. The low state 
and dominance of the companion over the spectra makes the 
classification of the system in terms of CV subclasses dif- 
ficult. An intermediate polar scenario, consistent with the 
presence of little disc, seems to be ruled out by the lack 
of typical lines from magnetic CVs in the spectra (such as 
prominent Hell, IWilliamsll 19891 ISchwarz et al.||2004l), but 



even those lines could be absent if the accretion rate is very 
low. The complex Ha line resembles nova-like systems such 
as BB Doradus ( Schmidtobreick et al.|20"l2 l which do show 
composite emission features that are not directly from ei- 
ther the secondary star or a remnant disc. Although the 
dynamics of the line suggests an origin near the tip of the 
companion star, in fact, the observed FWHM is too broad 
(~ 400 kms -1 ) to be ascribed to photospheric activity or ir- 
radiation of the secondary surface only (although saturation 
effect might play an important role). CX93 also shows some 
similarity with the VY Sculptoris subclass of nova-like CVs 
(NL), which do show extended low accretion states when 
there can be hardly any evidence of ongoing accretion (see, 
e.g. the low state of MV Lyrae in |Linnell et al.|2005[ ). A dwarf 
nova scenario also seems viable, with the source ongoing a 
small outburst at the time of the Mosaic-2 observations, but 
with the low accretion rate preventing clear accretion fea- 
tures from showing the classification remains tentative. 
CX93 also shows that a little level of variability in a quies- 
cence CV can produce spectral features resembling a quies- 
cent BH or NS. Plotting the EW of the Balmer lines detected 



in the VIMOS spectra against the X-ray to optical flux ratio 
(lag(F x /F v ) ~ -1, following [Patterson fc Raymond|p85 l 
provides outlying values for a CV that are suggestive of a 
more compact accretor. The values are consistent with a 
normal CV if we consider the IMACS spectra instead. More- 
over, it is worth noticing that reference quantities such as 
the X-ray luminosity and Fx/Fy are not good indicators 
of the source type for deeply quiescent systems as they are 
for actively accreting ones. The typical X-ray luminosity 
for quiescent NS or BH XRBs (10 31 - 10 32 ergs' 1 ) is in fact 
consistent with that of a CV in a low state, and Fy is the 
same for all XRBs when the companion star dominates the 
optical light. In conclusion, in a quiescent system like CX93 
even the nature of the accretor can be mistaken without a 
complete dynamical study (see also Marsh et al. 1994). This 
is important to take into account for projects, such as the 
CBS, that aim to identify new NSs or BHs. 



5 CONCLUSION 

We performed a full dynamical study of the CBS source CX- 
OGBS J174444.7-260330 (CX93), finding that the source is 
a nearby (400 < d < 700 pc) long orbital period (> 5 hours) 
CV, in a state of low accretion rate. One episode of accretion 
at higher rate was possibly caught in one of our lightcurves, 
displaying asymmetric maxima that could be fitted with a 
disc spot. Alternative causes for the asymmetric maxima, 
such as the presence of starspots on the companion star, 
are difficult to realistically model. The secondary is likely 
of spectral type (close to) K5 V and dominates the spectra, 
which display a disc contribution of ~ 2%. The classifica- 
tion of the source among CV subclasses is difficult. A DN 
(in particular a VY Sculptoris) or NL system in a low state 
seems likely, although we can not exclude a magnetic WD 
accreting at a very low rate. The case of CX93 highlights 
that CVs in a low accretion state can be hard to distinguish 
from quiescent XRBs with a NS or BH primary. 
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